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0.0  Abstract 

The  transformation  of  DOC  into  heterotrophic  biomass  is  an  important  ecological 
pathway.  This  paper  studies  the  theoretical  basis  for  variability  in  the  transformation 
resulting  from  aggregations  of  DOC  microzones  or  microclouds.  The  paper  computes 
measures  of  aggregation  using  the  theory  of  stochastic  geometry.  The  theory  enables 
computation  of  aggregation  in  terms  of  the  volume  fraction  occupied  by  microclouds  and 
computation  of  length  scales  such  as  covariance  functions  and  spherical  contact  distribu¬ 
tions.  Using  real  data  on  phytoplankton  cell  density  and  size,  and  conjectured  dimensions 
of  of  microzones,  we  were  able  to  compute  the  volume  fraction,  covariance  functions  and 
spherical  contact  distributions  for  assembledges  of  Baltic  Sea  phytoplankton.  By  compar¬ 
ing  microcloud  length  scales  with  molecular  diffusion,  turbulent  diffusion  and  uncorrelated 
velocity  length  scales  we  determined  that  variability  in  the  turbulent  kenetic  energy  dissi¬ 
pation  rate,  which  might  be  induced  by  global  and  basin  scale  wind  forcing  (for  example), 
could  influence  the  geometric  structure  of  carbon  microclouds.  This  demonstrates  the 
plausibility  of  linking  the  largest  scale  physical  forcing  directly  to  the  microscale  structure 
when  considering  variability  in  the  transformation  of  primary  to  secondary  production. 

1.0  Introduction 

An  important  component  of  oceanic  primary-secondary  production  transformation 
involves  the  absorption  of  dissolved  organic  carbon  (DOC)  by  bacteria.  There  are  sev¬ 
eral  sources  of  variability  associated  with  this  transformation.  One  source  of  variability 
involves  the  fraction  of  total  DOC  absorbed  as  exudates  at  relatively  high  DOC  concen¬ 
trations  very  close  to  the  photosynthesizing  phytoplankton  cells.  This  paper  developes  a 
theory  on  the  spatial  distribution  of  DOC  in  the  proximity  of  phytoplankton  cells.  The 
theory  generalizes  the  microzone  model  of  bioactive  DOC  distribution  near  metabolizing 
planktonic  organisms  (see  Mitchell  et  al.,  1985).  In  the  microzone  model,  exuded-carbon 
molecules  are  concentrated  in  a  spherical  zone  centered  on  each  phytoplankton  cell.  The 
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concentration  of  molecules  is  highest  at  the  cell  wall  and  then  declines  according  to  the  laws 
of  molecular  diffusion.  The  spherical  zones  are  thought  of  as  phy cospheres  or  microzones 
(The  idea  of  a  phycosphere  evidently  originated  in  Bell  and  Mitchell,  1972). 

However,  in  this  paper  we  postulate  that  the  microzones  tend  to  intersect  as  the 
numerical  density  of  the  phytoplankton  cells  increases.  A  particular  set  of  intersecting 
microzones  can  be  thought  of  as  a  “microcloud” .  The  ensemble  of  single  microzones  and 
microclouds  comprises  the  microcloud  structure-essentially  different  than  the  microzone 
structure-of  DOC.  Because  microclouds  are  larger  than  microzones,  they  are  subject  to 
different  modalities  and  intensity  of  physical  forcing.  The  following  sections  of  the  paper 
discuss  the  theory,  how  it  is  applied  and  a  specific  example. 


2.0  Theory 

Using  the  calculus  of  stochastic  geometry  (see  Stoyan  et  ah,  1987;  Cressie,  1993;  and 
Rothschild,  1992),  the  concentration  and  length-scale  properties  of  microclouds  can  be 
derived  as  a  special  case  of  the  general  Boolean  model.  The  model  enables  computation  of 
the  volume  fraction,  concentration  factor,  covariance  function  and  spherical  contact  dis¬ 
tribution  of  the  microclouds.  The  volume  fraction  and  the  concentration  factor  statistics 
are  measures  of  the  influence  of  microclouds  on  DOC  concentration  while  the  covariance 
function  and  the  spherical  contact  distribution  are  probability  distributions  reflecting  mi¬ 
crocloud  length  scales. 

These  statistics  and  distributions  are  developed  from  the  general  theory  by  assuming: 
(1)  that  the  phytoplankton  cells  are  distributed  in  5ft3  according  to  the  Poisson  distribution 
with  intensity  A  and  (2)  that  the  radii  of  the  microzones  are  distributed  log-normally  with 
mean,  £,  and  variance,  a2.  Accordingly,  the  probability  density  function,  /,  of  the  radii  is 
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where  exp  {•}  is  the  exponential  function  and 
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The  volume  fraction,  p ,  occupied  by  the  ensemble  of  microclouds  is, 


The  volume  fraction  described  in  equation  (2)  allows  for  the  possibility  that  the  mi¬ 
crozones  intersect  one  another.  The  intersection  of  microzones  can  serve  to  enhance  DOC 
concentrations  in  the  microcloud  over  those  in  a  microzone.  An  upper  bound  for  this  con¬ 
centration  effect  can  obtained  by  assuming  an  average  concentration  of  carbon  molecules 
within  the  sphere,  so  the  concentration  factor ,  Cf ,  is  the  ratio  of  the  volume  fraction 
without  overlap  to  the  actual  volume  fraction: 

^ttA  £3a3 

cf  =  ^— 

p 

The  concentration  factor  achieves  its  maximum  value  when  the  volume  fraction  without 
overlap  is  one.  This  yields  a  maximum  value  of  <7/max  =  (l  —  e-1)  1  w  1.58. 

The  length  scale  for  the  statistical  “size”  of  the  average  microcloud  can  be  derived 
from  the  stochastic-geometry  covariance  function.  The  covariance,  C(7Z),  is  the  probability 
that  both  the  head  and  tail  of  a  vector  of  length  1Z  are  located  within  the  volume  occupied 
by  the  microclouds.  Since  this  will  be  most  often  accomplished  by  having  the  vector  within 
a  single  microcloud,  the  covariance  can  be  thought  of  as  indicating  the  average  size  of  the 
microclouds.  Very  small  vectors  have  a  relatively  high  probability,  p,  of  being  contained 
within  the  average  microcloud,  while  large  vectors  haver  a  lower  probability,  p2,  of  being 
contained  within  the  microcloud.  The  covariance  function  is  given  by, 


C(1Z)  =  p2  +  (1  —  p)2  je3irA[(T’3'Ht1’  * ))  2,R-(r2H(r  * ))+ 16^  (H(r  *))]— lj 

where  (•)  denotes  an  expected  value  defined  by 
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A  convenient  rescaling  of  the  covariance  function  is  given  by  the  correlation  function : 


k(  %)  = 


cm  -  j_ 

p(i  -p) 


(5) 


Clearly,  k  (E  [0,1]  with  k  =  1  for  perfect  correlation  and  k  =  0  for  zero  correlation.  This 
makes  for  easy  comparisons  of  cases  with  different  volume  fractions,  as  well  as  correspond¬ 
ing  with  standard  conventions  for  the  notion  of  correlation. 

The  probability  that  a  point  selected  at  random  from  outside  the  volume  occupied  by 
the  microclouds  is  within  a  distance  H  of  a  microcloud  is  given  by  the  spherical  contact 
distribution ,  HS(1Z),  which,  under  our  assumptions,  is  given  by 


HS(H)  =  1  -  e-|7rA[3^(<T2+^2)+3'^2^+'^3] 


(6) 


3.0  Application 


In  order  to  place  equations  (l)-(6)  in  the  context  of  the  DOC  microcloud  problem, 
we  need  to  make  choices  for  A,  £,  a2  and  7 Z.  For  purposes  of  demonstration,  size-binned 
particle  counter  data  from  the  Baltic  sea  were  used  (Kahru  et  al .,  1991).  These  size 
binned,  or  size  frequency,  data  were  used  to  compute  estimates  of  the  numerical  cell 
density,  A  and  the  lognormal  parameters  £  and  a2.  The  intensity,  A,  ranges  between 
[lO3,  5  (lO3)]  (cells  •  cm-3).  These  spatially  averaged  concentrations  probably  represent 
much  higher  local  concentrations  because  of  the  known  non-random  distribution  of  phyto¬ 
plankton  cells. 

With  respect  to  the  mean  radius  of  the  microzone,  £,  direct  measurements  have  not 
been  made  (Mitchell  et  al. ,  1985,  have  conjectured  microzone  dimensions  and  experiments 
have  demonstrated  the  existance  of  (pH)  microzones,  Richardson  and  Stolzenbach,  1995). 
The  theory  of  diffusion  can,  however,  be  used  to  approximate  the  microzone  radius.  It  is 
well  known  that  the  steady  state  concentration  of  molecules  diffusing  away  from  a  cell  of 
radius  rcen  is  given  by 


c{r)  =  Coo  +  (Ccell 


(7) 


where  ccen  is  the  concentration  at  the  cell  wall  and  ca 0  is  the  far-field  (or  background) 
concentration.  The  microzone  associated  with  the  phytoplankton  cell  is  that  sphere  within 
which  the  concentration  exceeds  a  threshold  value,  Cthreshoid-  We  choose  the  threshold 
concentration  such  that  the  difference  (threshold  —  Cqo)  is  a  fixed  fraction,  6 ,  of  the  difference 
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(ccell  ~  Coo)-  The  mean  radius  is  chosen  to  be  that  radius  at  which  the  concentration  has 
fallen  to  cthreshold-  That  is,  find  £  such  that 

c(£)  =  Cthreshold  (8) 

=  Coo  T  8  (Ccell  Coo) 

Using  (7),  this  gives  the  mean  radius  as 

t  =  T-f  •  (9) 

In  this  paper  we  consider  8  =  0.1  which  gives  a  microzone  radius  10  times  that  of  the  phyto¬ 
plankton  cell.  Fitting  the  Kahru  et  al.  data  to  log-normal  distributions  via  a  \2  functional 
shows  the  mean  phytoplankton  cell  radius  to  fall  in  the  range  rceii  £  [4  (l0-4)  ,  8  (lO-4)] 
cm.  This  corresponds  to  a  mean  microzone  radius  in  the  range 

£  £  [4  (10“3)  ,  8  (10~3)]  cm  .  (10) 

The  magnitude  of  the  variance  in  the  microzone  radii  can  arise  either  from  the  vari¬ 
ability  of  the  phytoplankton  cell  radii  or  from  a  lack  synchrony  of  in  exudation  time.  For 
the  Kahru  et  al.  data,  the  variance  is  found  to  lie  in  the  range: 

<r2  €  [0.5£2  ,  0.9£2]  C  [8  (10“6)  ,  6  (10“5)]  cm2  .  (11) 

The  variance  and  the  mean,  of  course,  specify  the  shape  of  the  log-normal  distribution. 

The  correlation  function  and  the  spherical  contact  distribution  both  depend  upon  an 
additional  length  interval,  1Z  (either  a  correlation  distance  or  a  distance  to  the  nearest 
cloud).  The  range  of  7Z  is  determined  by  the  questions  being  asked.  We  wish  to  examine 
the  effect  of  physical  processes  (i.e.  turbulence)  on  the  microcloud  structure.  In  particular, 
we  consider  the  Kolmogorov  length  scale,  r],  and  the  turbulent  scalar  diffusion  length  scale, 
Tjs.  The  Kolmogorov  length  scale  is  the  smallest  length  scale  on  which  the  velocities  are 
uncorrelated,  while  the  turbulent  scalar  diffusion  length  scale  is  the  smallest  length  scale 
at  which  the  Kolmogorov  eddies  can  strain  a  tracer.  For  the  ocean,  typical  ranges  for 
the  Kolmogorov  and  turbulent  scalar  diffusion  length  scales  are  rj  £  [0.1,  1.0]  cm  and 
rjs  £  [l0~3,  0.1]  cm  (the  lowest  values  of  the  Kolmogorov  and  turbulent  scalar  diffusion 
length  scale  ranges  correspond  to  more  energetic  cases).  The  length  interval,  7 Z,  is  then 
chosen  to  contain  these  ranges: 

TZ  D  [10-3  ,  1]  cm  .  (12) 
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4.0  Results 

4.1  Estimates  of  Baltic  phytoplankton  statistical  geometry 
parameters. 

Fitting  the  Kahru  et  al.  data  to  the  log-normal  distribution  results  in  the  curves  of 
figure  1.  The  modal  microzone  radii  fall  in  the  range  of  [23,  33]  fxm.  The  curves  themselves 
span  the  range  of  those  with  small  variability  in  radii  to  those  with  larger  variability.  As 
a  general  trend,  the  populations  with  larger  radial  variability  are  those  with  larger  modal 
radii. 

The  resulting  statistics  are  more  conveniently  displayed  in  terms  of  nondimensional 
variables.  Noting  that,  on  average,  the  volume  occupied  by  one  particle  is  A-1,  then  a 
sphere  with  this  volume  will  have  a  radius,  R^,  given  by 


The  mean,  £,  and  length  interval,  7Z,  are  nondimensionalized  by  dividing  by  R*.  The 
variance,  <r2,  is  nondimensionalized  by  dividing  by  R^. 

Figure  2  shows  the  volume  fraction  as  a  function  of  the  relative  mean  radius,  and 
2 

the  relative  variance,  Circles  are  used  to  indicate  the  locations  in  parameter  space  of 
the  Kahru  et  al.  data,  with  a  bounding  region  for  the  data  shaded.  The  vertices  of  the 
bounding  region  are  labelled  (a)-(d)  for  comparison  with  figures  3-5.  The  volume  fraction 
occupied  by  this  Baltic  data  ranges  from  0.1%  to  10%. 

Figure  3  displays  the  concentration  factor  as  a  function  of  the  relative  mean  radius, 
-A-,  and  the  relative  variance,  j^-.  The  circles  and  shaded  region  are  exactly  the  same 
as  in  figure  2.  The  concentration  effect  in  the  shaded  region  is  bounded  between  a  0.03% 
increase  and  a  3%  increase.  This  reflects  that  there  is  relatively  little  overlap  in  the  Baltic 
phytoplankton  population. 

Figure  4  displays  contours  of  the  correlation  function  as  a  function  of  the  relative 
length,  and  the  relative  mean  radius,  The  variance  is  held  equal  to  half  the  mean 
radius  squared  in  figure  4a,  and  equal  to  0.9  times  the  mean  radius  squared  in  figure  4b. 
Both  figures  show  a  maximum  in  the  correlation  as  a  function  of  the  relative  mean  radius. 
In  figure  4a,  the  maximum  occurs  at  about  ^  =  0.70,  while  in  figure  4b  the  maximum 
shifts  down  to  around  =  0.55. 

Figure  5  displays  contours  of  the  spherical  contact  distribution  as  a  function  of  the 
relative  separation,  and  the  relative  mean  radius,  The  variance  is  held  equal  to 
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half  the  mean  radius  squared  in  figure  5a,  and  equal  to  0.9  times  the  mean  radius  squared 
in  figure  5b.  For  small  relative  mean  radius,  <  1,  the  contours  in  both  figures  appear 
independent  of  the  relative  mean  radius.  In  fact  this  region  corresponds  to  mean  radii  so 
small  that  the  spherical  contact  distribution  essentially  “sees”  dimensionless  points.  As  the 
relative  mean  radius  increases,  the  contours  tilt  to  the  left,  making  it  easier  to  be  within 
a  given  distance  of  the  microclouds.  Although  it  is  hard  to  see,  the  contours  in  figure  5b, 
(cr2  =  0.9£2),  tilt  more  than  those  in  figure  5a,  (cr2  =  0.5£2).  As  a  particular  example,  in 
figure  5a  the  contour  Hs(7l)  =  0.9  intersects  the  top,  =  10,  at  a  relative  separation  of 
=  0.005,  while  in  figure  5b  the  intersection  occurs  at  =  0.004. 

4.2  Influence  of  physical  processes  on  Baltic  phytoplankton. 

The  length  scales  (the  correlation  function  and  the  spherical  contact  distribution)  of 
the  microzones  are  roughly  in  the  molecular  diffusion,  turbulent  diffusion  and  Kolmogorov 
uncorrelate  velocity  length  scales.  The  molecular  diffusion  length  scale  is  not  subject  to 
variability  in  the  turbulent  kinetic  energy  dissipation  rate,  e.  However,  both  the  turbulent 
diffusion  and  the  Kolmogorov  uncorrelated  velocity  length  scales  are  subject  to  variations 
in  e.  Oversimplistically,  it  is  well  known  that  e  in  the  upper  ocean  varies  with  the  velocity 
of  the  wind  field  etc.  This  means  that  variations  in  the  global  wind  field  affect  e,  and  so  it 
is  important  to  determine  whether  the  global  wind  field  potentially  affects  the  pathways 
by  which  DOC  exuded  by  phytoplankton  cells  is  converted  into  secondary  production. 

In  order  to  judge  the  plausibility  of  the  linkage  between  physical  forcing  and  the 
microcloud  length  scales  we  have  constructed  rectangle  in  figures  4  and  5  that  delimit  the 
molecular  diffusion,  turbulent  diffusion  and  turbulent  velocity  length  scales. 

In  figure  4,  we  can  see  that  microclouds  are  generally  subject  to  molecular  diffusion 
and  the  straining  of  the  fluid  that  results  from  turbulent  diffusion.  On  the  other  hand  we 
can  also  see  that  the  microclouds  are  not  subject  to  forcing  by  uncorrelated  velocities.  This 
interpretation  is  distinct  from  the  viewpoint  that  only  microzones  exist  because  microzones 
are  only  subject  to  molecular  diffusion. 

In  figure  5,  we  see  that  the  interaction  between  points  in  the  void  space  and  material 
inside  the  microzone  is  influenced  by  both  turbulent  diffusion  and  uncorrelated  veloci¬ 
ties,  hence  the  interaction  is  affected  by  variability  in  e.  The  significance  of  the  effect 
of  turbulent  diffusion  might  involve  the  intensity  of  microcloud  blending  in  a  dynamic 
setting;  while  the  significance  of  the  effect  of  uncorrelated  velocities  is  that  there  would 
be  increases  in  encounter  rates  between  particles  not  in  a  cloud  and  particles  in  a  cloud 
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(see  Rothschilde  and  Osborn,  1988).  This  would  mean  that  encounters  between  secondary 
producers,  especially  microzooplankton,  and  bacteria  feeding  on  DOC  in  or  near  a  cloud 
could  be  enhanced  or,  more  importantly,  vary  with,  for  example,  the  velocity  of  the  wind 
field. 

5.0  Discussion 

5.1  The  role  of  microclouds  in  regulating  production. 

The  microcloud  landscape  is  an  important  source  of  variability  in  regulating  the  cou¬ 
pling  between  primary  and  secondary  production.  The  notion  of  the  microcloud  suggests 
that  the  DOC  available  to  bacteria  may  be  higher  than  might  be  assumed  taking  account 
of  only  some  background  level  of  carbon.  In  addition  the  rich  variability  in  spatial  dis¬ 
tribution  implies  that  spatial  distribution  is  an  important  source  of  the  variability  in  the 
primary  to  secondary  production  via  the  bacterial  pathway. 

The  nature  of  the  spatial  coupling  might  be  best  thought  of  in  terms  of  an  initial 
random  distribution  of  bacteria.  Such  a  distribution  might  exist  if  carbon  molecules  were 
distributed  at  random.  However,  the  microzone  or  microcloud  notion  ensures  a  non-random 
distribution  of  carbon  molecules.  As  a  result  of  possible  tactic  responses  and  accelerated 
growth  rates  in  the  regions  of  relatively  high  carbon  molecule  concentration  it  can  be 
expected  that  the  densities  of  bacteria  in  the  microclouds  would  be  higher  than  outside 
the  microclouds.  Thus  high  volume  fractions  would  be  conducive  to  a  higher  transfer  rate 
between  primary  production  and  secondary  production. 

An  alternative  view  of  this  coupling  is  presented  in  the  work  of  Bowen  et  al.  (1993). 
In  their  simulations  a  single  phytoplankton  cell  was  centered  in  the  computational  volume 
and  1000  bacteria  were  initially  randomly  placed  around  it.  A  time  varying  shear  was  im¬ 
posed  on  the  fluid,  representing  the  shear  induced  by  turbulence.  Their  results  show  that, 
although  an  individual  bacterium  could  not  remain  in  the  neighborhood  of  a  phytoplank¬ 
ton  cell,  their  transient  proximity  was  sufficient  to  give  significantly  enhanced  exposures 
to  exudates.  This  could  be  a  very  important  mechanism  for  the  Baltic  sea  populations, 
were  the  contact  distribution  shows  that  it  is  very  easy  to  be  close  to  a  microcloud. 

It  would  follow  that  higher  concentration  of  bacteria  in  microzones  would  induce  higher 
concentrations  of  microzooplankton  and  protists  that  feed  upon  bacteria.  In  strengthing 
the  coupling  to  an  even  greater  degree,  the  increased  densities  of  the  microzooplankton 
would  also  result  from  tactic  respones  and  increased  growth  rates  which  are  known  to  be 
exceptionally  high  with  abundant  food  resources  (e.g.  Fenchel  1988). 
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The  aggregation  of  a  variety  of  small  heterotrophs  within  the  microcloud  of  carbon 
molecules  suggests  than  the  carbon  microcloud  is  juxtaposed  with  an  ammonium  micro¬ 
cloud  produced  by  the  heterotrophs,  closing  the  feedback  loop  and  increasing  to  yet  a 
greater  degree  the  strength  of  the  coupling.  The  synergy  of  this  coupling  would  have  the 
effect  of  enhancing  phytoplankton  growth  through  the  increased  availability  of  nutrients 

This  clearly  evokes  a  very  complex  structure.  On  one  hand,  from  the  phytoplankton’s 
point  of  view,  the  carbon  gradients  both  overlap  and  decrease  with  distance  from  the  cell 
wall.  On  the  other  hand,  because  of  diffusion,  the  ammonium  molecules  increase  in  con¬ 
centration  with  distance  from  the  phytoplankton  cell  wall.  Thus,  from  the  phytoplankton’s 
point  of  view,  the  environment  consists  of  both  a  positive  and  a  negative  microcloud.  The 
reverse  is,  of  course,  true  from  the  heterotroph  perspective. 

5.2  Negative  microclouds. 

The  negative  microcloud  is  particularly  important  because  its  existance  impacts  upon 
the  issue  of  diffusion  limited  nutrient  availability  for  phytoplankton  cells.  Munk  and  Riley 
(1952)  and  Pasciak  and  Gavis  (1974)  have  studied  this  problem  in  the  context  of  molecular 
diffusion.  Lazier  and  Mann  (1989)  observed  that  molecular  diffusion  is  influenced  at  very 
small  scales  by  turbulent  flow  and  hence  follows  turbulent  diffusion  at  very  small  scales.  To 
some  extent,  variable  concentration  may  be  “smoothed”  by  strain  resulting  from  turbulent 
flow.  At  any  rate,  all  of  the  above  results  need  to  be  modified  by  the  existance  of  negative 
microclouds.  In  effect,  the  negative  microcloud  generates  a  structure  of  rather  complex 
gradients  and  hence  complex  diffusional  flux.  In  fact,  it  is  possible  for  nutrient  molecules 
to  be  “attracted”  to  lacunae  inside  the  cloud  as  well  as  to  cells  actively  absorbing  nutrient 
molecules.  The  net  over- all  effect,  however,  would  be  to  generally  decrease  the  nutrient 
delivery  to  the  phytoplankton  cells  within  the  microcloud,  in  direct  competition  with  the 
nutrient  enhancement  described  in  the  preceeding  section.  The  net  result  of  these  two 
opposing  effects  is  an  important  issue  that  is  beyond  the  scope  of  this  paper. 

5.3  Future  directions. 

Thus  we  have  established  the  plausibility  of  the  existence  of  a  microcloud  structure 
mediating  the  transformation  of  exuded  DOC  into  secondary  production.  Perhaps  more 
importantly,  we  have  shown  how  this  transformation  can  be  modulated  by  physical  forcing 
associated  with  the  turbulent  kinetic  energy  dissipation  rate.  This  is  of  course  remarkable 
because  it  links  very  large  scale  processes  such  as  the  wind  field  directly  with  the  smallest 
scale  processes  that  relate  to  diffusion  and  turbulent  flow. 
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Future  research  would  include  the  development  of  laboratory  and  field  experiments 
to  test  the  theory.  Additional  theoretical  developments  would  include  understanding  how 
the  geometry  of  clouds  is  affected  by  the  dynamics  of  the  flow  field  and  the  properties  of 
diffusion. 

5.4  Acknowledgements. 
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ire  1  Probability  density  curves.  The  data  of  Kahru  et  al.  (1991)  fit  to  the  log 
.istribution.  Each  curve  refers  to  a  10  day  period  in  the  data. 
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Figure  2  Volume  fraction  contours.  The  volume  fraction  occupied  by  the  ensemble 
microzone.  The  circles  indicate  the  position  in  parameter  space  for  the  Kahru  et  ah  (1991) 
data.  The  shaded  region  is  the  “reference”  area  suggested  by  the  Kahru  et  ah  data. 


Figure  3  Concentration  factor  contours.  An  upper  bound  for  the  concentration  effect 
due  to  the  overlapping  of  the  microzones.  The  circles  indicate  the  position  in  parameter 
space  for  the  Kahru  et  al.  (1991)  data.  The  light  gray  shaded  region  is  the  “reference”  area 
suggested  by  the  Kahru  et  al.  data.  The  dark  gray  shaded  regions  in  the  upper  corners 
are  parameter  regimes  in  which  the  microzones  must  overlap. 
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Figure  4  Correlation  contours  for  (a)  a2  =  0.5£2  and  (b)  cr2  =  0.9£2.  Dotted  lines 
delimit  ranges  for  phytoplankton  microzone  mean  radii  of  the  Kahru  et  al.  (1991)  data, 
corresponding  to  points  (a)-(d)  in  figure  2. 
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Figure  5  Spherical  contact  distribution  contours  for  (a)  <r2  =  0.5£2  and  (b)  tr2  =  0.9£2. 
Dotted  lines  delimit  ranges  for  phytoplankton  microzone  mean  radii  of  the  Kahru  et  al. 
(1991)  data,  corresponding  to  points  (a)-(d)  in  figure  2. 


